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Abstract 

Binary radio pulsar system J0737— 3039 provides an exceptional opportunity to 
study innermost structure of pulsar magnetospheres due to very tight orbit, favorable 
directions of pulsars' rotation and magnetic axes and extremely fortuitous orientation 
of the orbit. In this system the millisecond pulsar A is eclipsed once per orbit. 
During eclipse a clear modulation at the 2.77 s period of pulsar B is seen, pointing 
unambiguously to magnetospheric origin of eclipses. A simple geometric model, 
based on the idea that the radio pulses are attenuated by synchrotron absorption on 
the closed magnetic field lines of pulsar B, can successfully reproduces the eclipse 
light curves down to intricate details. This detailed agreement confirms the dipolar 
structure of the neutron star's magnetic field. The model gives clear predictions for 
temporal evolution of eclipse profile due to geodetic precession of pulsar B. 

In addition, pulsar B shows orbital modulations of intensity, being especially bright 
at two short orbital phases. We showed that these modulations are due to distortion 
of pulsar B magnetosphere by pulsar A wind which produces orbital phase-dependent 
changes of the direction along which radio waves are emitted. Thus, pulsar B is 
intrinsically bright at all times but its radiation beam misses the Earth at most 
orbital phases. 
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1 Amazing and lucky: double 
pulsar PSR J0737-3039A/B 



The double pulsar system PSR 
J0737— 3039A/B contains a recy- 
cled 22.7 ms pulsar (A) in a 2.4 hr 
orbit around a 2.77 s pulsar (B) 
20031 : iLvne et al 



( Bur gav et al. 



edge on, with the line of sight making 
an angle ~ .5 ° with t he orbital plane 



( Col e et all 12004 iRansom et al. 



20041 ). This leads to a number of ex- 



2004). The system is viewed nearly 



ceptional properties of the system. 
Pulsar A is eclipsed once per orbit, 
for a duration of ~ 30 s centered 
around superior conjunction (when 
pulsar B is between the observer and 
pulsar A). The width of the eclipse is 
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only a wea k function of t he ob serving 
frequency ([Kaspi et all 120041 ) . Most 
surprisingly, during eclipse the pulsar 
A radio flux is modulated by the ro- 
tation of pulsar B: there are narrow, 
transparent windows in which the 
flux from pulsar A rises nearly to the 



unab sorbed level ([McLaughlin et al 



|2004| ). These spikes in the radio flux 
are tied to the rotational phase of 
pulsar B, and provide key constraints 
on the geometry of the absorbing 
plasma. At ingress transparent spikes 
appear at half rotational period of 
B, then change their modulation to 
rotational period of B at the middle 
of eclipse and disappear completely 
right before the egress (Fig 1). On 
average the eclipse is longer when 
the magnetic axis of pulsar B is ap- 
proximately aligned with the line of 
sight (assuming that radio pulses are 
generated near the magnetic axis of 
B). 

The physical width of the region 
which causes this periodic modula- 
tion is comparable to, or somewhat 
smaller, than the estimated size of 
the magnetosphere of pulsar B. Com- 
bined with the rotational modula- 
tion, this provides a strong hint that 
the absorption is occurring within 
the magnetosphere of pulsar B. This 
allows us to probe directly the struc- 
ture of pulsar B magnetosphere. 



presence of a dipolar magnetic field 
around pulsar B. 



2 Magnetosphere of pulsar B 



At distances much larger than neu- 
tron star radius but much smaller 
than light cylinder radius, Rlc ~ 
c/Q (Q is rotation frequency of a 
pulsar, magnetospheres of isolated 
pulsars are nearly dipolar, contain- 
ing two types of magnetic field lines: 
those that close inside Rlc and 
those that cross the light cylinder 
and thus remain o p en to infinity 



1969). Radio 



f)Goldreich fc Julian , 
emission is, presumably, generated 
on the open field lines near the neu- 
tron star magnetic polar field line. 

In case of PSR J0737-3039B this 
picture is modified by strong interac- 
tion of pulsar B with relativistic wind 
flowing outward from pulsar A. Wind 
of A strongly distorts B magneto- 
sphere, blowing off a large fraction of 
it, so that the size of magnetosphere 
of B facing A is much smaller than 
light cylinder and is now determined 
by the pressure balance between the 
wind of pulsar A and the magnetic 



press ure of pulsar B (e.g. iLvutikov , 
120041 ) . Observationally, this is con- 
firmed by the fact that the half- width 
of the eclipse, ~ 1.6 x 10 9 cm, is 
We constru cted a detailed model much smaller than the pul sar B light 

cylinder Rlc ~ 1-3 x 10 10 cm. 

As a results of interaction with A 
wind, spindown torque on pulsar 
B is modified if compared with the 
vacuum case. This has important 
implications for pulsar B magnetic 



of ec lipses fjLvutikov Thompson 
2005b which shows that the light 



curve is consistent, in considerable 
detail, with synchrotron absorp- 
tion in magnetosphere of pulsar B. 
The agreement is detailed enough 
to constitute direct evidence for the 
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field and spindown age, which can 
be estimated self-consistently only 
by modeling the interaction between 
wind and magnetosphere (conven- 
tional vacuum spindown formulas 
are incorrect in th is case!). Esti- 
mates of the torqu e (|Lvutikovl . 12004 ; 
Arons et all 120041 ) show that break- 
ing index of pulsar B is likely to be 
close to 1 and may even approach 
(instead of vacuum dipole case of 3). 
Surface magnetic field and magneto- 
spheric radius can be estimated as 
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where Qa,b and Qa,b are mea- 
sured spin fre quencies and f re quenc y 



derivatives ([Burgav et al. 



^quency 
2003), 



Ia,b are moments of inertia nor- 
malized to 10 45 g cm 2 and Dab is 
distance between pulsars. 

As long as parameters of pul- 
sar A remain approximately con- 
stant, over time scale of 250 Myrs, 
pulsar B spins down exponen- 
tially, Q B oc exp{— t/r} where 
r = Vt B /tl B « 100 Myrs (note that 
this is two times larger than its char- 
acteristic age). 

While details of magnetosphere-wind 
interaction are bound to be compli- 
cated, a reasonable guess, on which 



the above estimates are based, is 
that open-field current density re- 
mains similar to the one in isolated 
pulsars, j ~ ricjec {n GJ = QB/2irec 
is Goldreich- Julian density), while 
the size of open field lines is in- 
creased due to smaller size of B 
magnetosphere. In fact, field drag- 
ging at the magnetospheric bound- 
ary may increase typical current 
density deep inside magnetospheres. 
This will lead to increased torque 
and decreased estimates of surface 
magnetic field and magnetospheric 
radius. Still, we can estimate the 
minimum magnetospheric radius 
by assuming that at the boundary 
toroidal magnetic field does not ex- 
ceed poloidal magnetic field. This 
gives a minimum estimate of the 
magnetospheric radius (correspond- 
ing to spindown index 0): i? mag min = 

10 9 cm. Note, that this is still some- 
what larger than the size of eclipses. 



2. 1 Model of eclipses 



Similar to the case of isolated pul- 
sars, we assume that magnetic fields 
lines of pulsar B are of two types: 
open and closed ones. In addition, we 
assume that closed magnetic fields 
lines extend not up to the light cylin- 
der but up to the magnetosphere ra- 
dius -R mag , while field lines extending 
further than -R mag remain open. We 
place the absorbing relativistically 
hot plasma within a set of closed 
dipolar field lines. It can be shown 
that plasma density is constant along 
each field line. We further assume 
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that density and temperature do not 
vary between field lines. Both of 
these assumptions have only minor 
effect on eclipse profile. 

We calculate the synchrotron optical 
depth over a large number of lines of 
sight, taking into account both the 
three-dimensional structure and ro- 
tation of B magnetosphere (Fig. 2). 
The calculation of the eclipse light 
curve requires a choice of several pa- 
rameters: direction of rotational axis 
(angles 9q, 4>q), angle between rota- 
tional and magnetic axis Xi impact 
parameter zq, electron temperature 
T e , density multiplicity A mag (ratio of 
density to Goldreich- Julian density) 
and outer radii -R a bs+ 

Results of the simulations are com- 
pared with data on in Fig. (1) The 
model reproduces many fine details 
of the eclipses. The model can re- 
produce asymmetric form of eclipse 
with the ingress shallower and longer 
than the egress. It explains the mod- 
ulation of that is observed at the 
first and second harmonics of the 
spin frequency of pulsar B, and the 
deepening of the eclipse after supe- 
rior conjunction. The average eclipse 
duration is almost independent of 
frequency when multiplicity is suf- 
ficiently large. The onset and ter- 
mination of the eclipse are deter- 
mined mostly by the physical bound- 
ary of the absorbing plasma and 
not by microphysics of absorption 
process. This results in nearly fre- 
quency independent eclipses. Eclipse 
is broader when the magnetic mo- 
ment of pulsar B is pointing closest 
to the observer, ju st as is observed 



are modest deviations between the 
model and the data near the edges 
of the eclipse that could be used to 
probe the distortion of magnetic field 
lines from a true dipole. 

The modulation of the radio flux 
during the eclipse is due to the fact 
that - at some rotational phases 
of pulsar B - the line of sight only 
passes through open magnetic field 
lines where absorption is assumed to 
be negligible. One of the main suc- 
cesses of the model is its ability to 
reproduce both the single and dou- 
ble periodicities of these transparent 
windows, at appropriate places in 
the eclipse. This requires that /2g 
be approximately - but not quite - 
orthogonal to fi^. 

Based on a sample of many eclipse 
light curves our best fit parameters 
are: z ~ —7.5 x 10 8 cm, 9q ~ 60°, 
<p n ~ -90°, x - 75°. The orbital in- 
clination is therefor e predicted to be 
close to 90.55° , (cf. " 



Cole et all 12004: 



Ransom et all 120041 ) . The size of the 
eclipsing region is -R a bs+ = 1.5 x 10 9 ; 
the implied multiplicity is large, 



A 



n 



mai] 



/n GJ 



10 . Lorentz 



fjMcLaughlin et all 120041 ) . There 



factor of suspended plasma was fixed 
at 7 = 10 (there is a degeneracy 
between Lorentz factor and density) . 

Presence of high multiplicity, rela- 
tivistically hot plasma on closed field 
lines of pulsar B is somewhat surpris- 
ing, but not unreasonable. Dense, 
relativistically hot plasma can be 
effectively stored in the outer mag- 
netosphere, where cyclotron cooling 
is slow. The gradual loss of particles 
inward through the cooling radius, 
occurring on time scale of millions of 
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pulsar B periods, can be easily com- 
pensated by a relatively weak upward 
flux driven by a fluctuating compo- 
nent of the current. For example, if 
suspended material is resupplied at a 
rate of one Goldreich-Julina density 
per B period and particle residence 
time is million periods, equilibrium 
density will be as high as 10 6 ncj- 
The trapped particles are also heated 
to relativistic energies by the damp- 
i ng of magnetospheric turb ulence 
( Lvutikov fc Thompson . 2005f ). 



There is only one slight deviation 
from simple estimates: our eclipse 
calculations show that the optical 
depth of the absorbing plasma un- 
dergoes a sharp drop at a distance 
-Rabs+ — 1-5 x 10 9 cm from pulsar B. 
This is about 2.5 times smaller than 
the expected radius of the magne- 
topause, -R mag ~ 4 x 10 10 cm. Most 
likely this is due to the loss of plasma 
from the outermost closed field lines 
due to reconnection and/or gradient 
drift. 



2.2 Predictions of the model 



There is a number of predictions of 
the model that should be tested in 
the coming years. 

(1) The spin tts of pulsar B is ex- 
pected to undergo geodetic pre- 
cession on a ~ 75 yea r timescale 
( Burgavet all . l2003h . We have 
provided predictions for how 
the eclipse light curve will vary 
as a result. In particular, the 
orbital phase at which the ra- 
dio flux reaches a minimum will 



shift back toward superior con- 
junction. Since (fin is not well 
constrained, a time for eclipse to 
become symmetric is between 
~ 12 years (if \(fi n + 90° | ~ 30° 
and Qb is drifting away from 
the plane of the sky) and ~ 25 
years (if (fin is drifting toward 
the plane of the sky) . 

(2) High time-resolution observa- 
tions are a sensitive probe of 
the distribution of plasma prop- 
erties (density and tempera- 
ture) on the closed field lines. If 
plasma is depleted from the out- 
ermost field lines, at high tem- 
poral resolution the flux should 
return to unity. On the other 
hand, if the absorbing plasma 
does not have a sharp trunca- 
tion in radius, then the flux will 
not return to unity in all of the 
transparent windows. 

(3) The eclipses must regain a 
strong frequency dependence 
at sufficiently high frequen- 
cies. The critical frequency 
above which significant trans- 
mission occurs can be used to 
place tight constraints on the 
plasma density. The electron cy- 
clotron frequency at a distance 
r ~ zq ~ 7.5 x 10 8 cm from pul- 
sar B is estimated to be Vs,e ~ 3 
GHz. One therefore expects the 
eclipses to develop a significant 
frequency dependence at higher 
frequencies. 

(4) There number of definite 
predictions for the polarization 
of the transmitted radiatio n 
( Lvutikov fc Thompson . 20051 ). 
But since pulsar A emission is 
strongly polarized one needs to 
know absolute position of the 
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direction of linear polarization. 



3 Orbital modulation of B due 
to distortion of B magneto- 
sphere by A wind 



Next we turn to orbital modulations 
of intensity of J0737-3039B. Emis- 
sion of pulsar B is strongly dependent 
on orbital phase, being especially 
bright at two wi ndows, each lastin g 



for about 30° fjLvne et all 120041) 



Pulse profiles have different shapes 
i n the two wi ndows. We showed 
f)Lvutikovl . l2005h ihat orbital bright- 
ing of B is not intrinsic but is due 
to distortions of B magnetosphere 
by pulsar A wind, so that radiating 
beam of B mostly misses the ob- 
server, except at two orbital phases. 

To prove this point, we assume that 
intrinsically bright pulsar B radio 
emission is generated near the polar 
field line inside a region with half 
opening angle of ~ 2 degrees, while 
the polar field line itself is "pushed 
around" by the pressure of A wind. 
To parametrize the effects of A wind 
on direction of radio emission of B 
we employ a method of distortion 
transfo r matio n of Euler potentials 



(Stern], Il994f l . developed for mod- 



eling of the Earth magnetic field 
influenced by solar wind. Assuming 
that radio emission of B is generated 
in physically small region of B mag- 
netosphere (much smaller than its 
size), direction of radio emission of 
B can be parametr ized by a distor- 
tion coefficient C ( Lvutikov . 20051 ) 
(so that distortion is proportional to 
1 — C, for pure dipole C = 1). 



To fit the data we use the results 
of modeling of pulsar A eclipse 
( Lvutikov fc Thompson , 2Q05h . Search- 
ing through parameter space we 
were able to obtains a satisfactory 
fit, producing emission at two or- 
bital phases nearly coincident with 
observed ones, Fig. 3. 



The strength of the model is that 
it reproduces (or predicted) fairly 
well a number of properties of pul- 
sar B radio emission. First, since at 
different orbital phases the line of 
sight makes different cuts through 
emission region, the model naturally 
explains orbital phase- dependent av- 
eraged emission profile. In fact, this 
can be used to construct a map of the 
emission region, something which 
was not possible for isolated pulsars. 
Secondly, since the emitting beam 
of B is pointing in somewhat dif- 
ferent directions at different orbital 
phases, time of arrival of pulses of 
B should experience regular orbital- 
dependent changes by as much as 
~ (2°/360°)2.8sec ~ 10 ms. Such 
drift in t imes of arriva l has been 
observed: iRansom et al. (2004) re- 
ported a systematic change in arrival 
times of B pulses by 10-20 ms. 



Finally, since emission beam of B is 
fairly narrow, small changes in the 
orientation of rotation axis of B may 
induce large apparent changes in the 
profile. One possibility for the change 
in the direction of the spin is geode- 
tic precession of B, which should hap- 
pen on a relatively s hort t ime scale ~ 



70 yrs ([Lvne et all 120041 ) . From our 



modeling we find that changes of <J)q 
by as little as ~ 1° strongly affect ob- 
served B profile. Thus, we expect that 
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profile of B may change on a times 
scale of less than a year. 1 This pre- 
di ction has been recen tly confirmed 
by Burgav et al. ( 200ot , who see evo- 
lution of profile of PSR J0737-3039B 
in general agreement with our model. 

Another possible implication of the 
model concerns emission heights in 
pulsar B. Data are best fitted by the 
model with the stretching coefficient 
C = 0.7, so relative deformation of 
the magnetosphere at the emission 
radius is ~ 30%. This is a fairly large 
distortion, which favors large emis- 
sion altitudes, R em > 10 8 cm. Large 
emission altitudes in isolated pulsars 
have been previously suggested by 
Lvutikovet aD(|l998l) . 



structure of magnetosphere may de- 
pend on the direction of the wind 
magnetic field. We plan to address 
these issues in a subsequent paper. 



The main implication of the model is 
that B is always intrinsically bright. 
Its peak luminosity of ~ 3 mJy at 
820MHz is, in fact, typical for iso- 
lated pulsars with similar properties. 
There is a number of ways the model 
can be improved. First, a better 
model of dipole distortion is needed. 
This can be achieved by using well 
developed models of the Earth mag- 
netosphere: since the size of pulsar B 
magnetosphere is much smaller than 
light cylinder, light travel effects are 
probably unimportant. Second, a 
non-trivial geometry of the emission 
region, as opposed to a circular used 
in the present study, may increase 
the quality of the fit. Also, if recon- 
nection between wind and magneto- 
spheric field lines is important, the 



1 According to the model, we are not 
likely to lose pulsar B in the coming 
year(s), yet the model is not sufficiently 
detailed to guarantee it. 



4 Conclusion 



The success of the models is some- 
what surprising, given the complex- 
ity of data and the fact that we had 
to fit many parameters. The reasons 
for it is that, to the first order, mag- 
netic field of B is well approximated 
by dipolar structure. Thus, the mod- 
els provide a strong test of the long- 
standing assumption that neutron 
stars are surrounded by corotating, 
dipolar magnetic fields. This is a 
valuable confirmation of a funda- 
mental assumption made in models 
of pulsar electrodynamics. In addi- 
tion, our results are consistent with 
models that place the source of the 
radio emission close to the magnetic 
axis. Forthcoming improvement of 
the model should be able to quan- 
tify deviations of the magnetosphere 
from dipolar shape and probe how a 
pulsar magnetosphere interacts with 
an external wind. 

Our eclipse model is consistent with 
a large natal kick of pulsar B that 
changes the orientation of the orbital 
plane and disrupt any pre-existing 
alignment between orbit and the spin 
of the progenitor star. 
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Orbital phase (degrees) 

Fig. 1. Comparison of a simulated 
eclipse profile with 800 MHz data 
( McLaughlin et all 12004^ 1 . The model 
fits the data best in the middle of the 
eclipse. 
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2. View of the magnetosphere at different rotational phases separated by tt/4. 
a full movie of the eclipse see http:/ /www. physics. mcgill.ca/~lyutikov/movie.gif 
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indicate orbital phases where B emission is strongest, dark regions indicate the 
location of emission in our best fit model. 
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